Astronomy & Astrophysics manuscript no. 14685 
May 24, 2010 



©ESO2010 



Letter to the Editor 



o 

(N 



o 

(N 

o 
u 

^' 

Oh 
I • 

o- 



First detection of tlie Sunyaev ZeI'dovich effect increment at 

A < 650jL/m* 

M. Zemcov[22], M. ReM, T. D. RawlJ^l, J. J. Boc]i221, E. EgamP, B. AltierJ^l, A. W. Blaiiff, R BoonJSEl^ 

C. R. Bridg^R ClemenP, R CombeJ^l, c. D. DowelffEl, M. Dessauges-ZavadskjlS], D. FaddtP, o. Ilber^, 

R. J. IvisoiffSHI, M. JauzacP, J.-P. KneibftD. LutJEl, R. Pell(J2, M. J. Pereiri^l, P G. Perez-Gonzale^l^, 

J. RicharcfH], g. H. RiekJ^l, G. RodighiercffS], d. SchaereJSEl, G. P SmitlJIS, I. Valtchano\Sl, G. L. WaltlP, 

P van der Werli^l, and M. W. WerneP 

(Affiliations can be found after the references) 

Submitted April 1, 2010, accepted May 10, 2010 

ABSTRACT 

The Sunyaev-Zel'dovich (SZ) effect is a spectral distortion of the cosmic microwave background as observed through the hot plasma in galaxy 
clusters. This distortion is a decrement in the CMB intensity for A > 13 mm, an increment at shorter wavelengths, and small again hy A ~ 250 yum. 
As part of the Herschel Lensing Survey (HLS) we have mapped 1E0657-56 (the Bullet cluster) with SPIRE with bands centered at 250, 350 
and 500 yum and have detected the SZ effect at the two longest wavelengths. The measured SZ effect increment central intensities are A/o = 
0.097 ± 0.019 MJy sr"' at 350 /jm and A/q = 0.268 ± 0.031 MJy sr"' at 500 /jm, consistent with the SZ effect spectrum derived from previous 
measurements at 2 mm. No other diffuse emission is detected. The presence of the finite temperature SZ effect correction is preferred by the SPIRE 
data at a significance of 2. Icr, opening the possibility that the relativistic SZ effect correction can be constrained by SPIRE in a sample of clusters. 
The results presented here have important ramifications for both sub-mm measurements of galaxy clusters and blank field surveys with SPIRE. 

Key words, cosmic background radiation - Galaxies: clusters: individual: IE0657-56 
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1. Introduction 

The Sunyaev-Zel'dovich (SZ) effect is a distortion of the 
spectral shape of the cosmic microwave background (CMB) 
due to inverse Compton scattering in the ubiquitous, hot 

I Je ~ 10^ K) intrac luster medium (ICM) of galaxy clusters 
Sunyaev & Zeldovichlll972 ). The canonical thermal SZ spec- 
trum is a decrement in the brightness of the CMB as measured 
through galaxy clusters at mm wavelengths and an increment at 
sub-mm wavelengths which passes though a null ?A A oi 1.3 mm. 
To correctly describe the SZ spectral distortion when relativis- 
tic electrons are present or the cluster is moving with respect to 
the CMB additional correction terms, usually termed "relativis- 
tic" (or "finite temperature") coiTections and "kinetic" SZ effect, 
are required. Measurement of these coiTections is only possi- 
ble using observations at multiple wavelengths, and is expedited 
by measurement at wavelengths where the expected modifica- 
tions to the thermal SZ effect are largest. In the case of the finite 
temperature corrections, the largest changes expected in the SZ 
increment are at wavelengths shorter than 1 /vm. 

Decrements in emission are rare astrophysically and can 
be ascribed to the SZ effect with little ambiguity; this has 



* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. It is open for proposals for observing 
time from the worldwide astronomical community. Data presented in 
this paper were analyzed using "The Herschel Interactive Processing 
Environment (HIPE)," a joint development by the Herschel Science 
Ground Segment Consortium, consisting of ESA, the NASA Herschel 
Science Center, and the HIFl, PACS and SPIRE consortia. 



lead to measurements of the SZ effect ai A > 2 mm be- 
coming almost routine. In comparison, measurements of the 
SZ effect increment are complicated by the presence of the 
dusty, high redshift galaxies which constitute the sub-mm cos- 
mic background. Additionally, the individual sources compris- 
ing the sub-mm background are gravitationally lensed by galaxy 
clusters, the effect of which is to preferentially coiTelate in- 
creases in sub-mm emission with clustering. This coiTelation 
makes unambiguous detection of the SZ increment difficult. 
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Komatsu et alj |l999l IZemcov et all I2003L , , . 

,Nord et alj|2009h . Moreover, the presence of the sub-mm back- 
ground may contaminate measurements of the SZ effect fo r 
A> \ mm in less massive galaxy clusters dAghanim et al.ll2005h . 
A better understanding of the sub-mm emission associated with 
galaxy clusters is required. 

Heretofore, systematic far infrared (FIR) surveys of many 
galaxy clusters to large radii have been technically challenging 
so a complete census of sub-mm emission f rom clusters has been 
difficult to obtain. T he advent of SPIRE ( Griffin etal.1 12010) on 
Herschel dPilbratt et al. 2010) has, for the first time, provided 
the capability to make deep maps of clusters to large angles on 
the sky and to use colour information to separate the different 
sources of sub-mm emission present in galaxy clust ers. In ad- 
dition to gravitationally lensed background sources dRex et alj 
2010), emission in clusters above the confused sub-mm back- 
groun d may also compri se emission from galaxies in the cluster 
itself dRawIe et alj|2010l) . as well as truly diffuse emission from 
the SZ effect and possibly even cold dust in the ICM. SPIRE's 
abiUty to separate sources of emission based both on spatial and 
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spectral information allows the demographics of the sub-mm 
emission to be measured. 

In this paper, we use deep SPIRE maps of the z - 0.3 Bullet 
cluster (1E0657-56) taken as part of the Herschel Lensing 
Survey (HLS, P.I. Egami) at 250, 350, and 500 ^/m with 18, 25 
and 36 arcsec resolution to measure the SZ effect and constrain 
other diffuse emission associated with the cluster. 

2. Analysis and results 

The HLS is a programme to observe 40 massive clusters in the 
range 0.1 < z < 1.0; during Herschel's science demonstra- 
tion phase 1E0657-56 was observed to full depth. The specific 
observation strategy, low le vel data reduc tion and map making 
processes are summarized in lEgami et alj (2010); here we begin 
with the calibrated flux and noise maps at 250, 350 and 500 //m 
(Fig. 2 of lEgamietani2010h . 

Measurement of dim, extended emission is complicated 
by the presence of background sub-mm galaxies the confu- 
sion noise from which, at SPIRE's resolution, i s measured to 
be ~ 6mJy in all three bands for blank fields (" Nguyen et alj 
I2OIOI) . Indeed, in these maps we measure an RMS noise of 
{6.0, 6.0, 7.5) mJy/beam at 250, 350, 500 yum, consistent with the 
combination of confusion noise and instrument noise expected 
from the integration time. Obviously, care must be taken to dis- 
entangle extended emission from confused sources. To do this, 
we employ the approach discussed below, the fundamental tenet 
of which is to account for and remove all 250 //m emission from 
the two longer wavelength channels. 

To generate a 250 yum source catalog, sources are detected in 
the 250 yum map using the StarFinder algorithm dPiolaiti et alj 
[2000). The detection threshold is set to 2.5cr and is run itera- 
tively using SPIRE's 250 /vm point response function (PRE) to 
account for blending; as the purpose here is not to generate a 
catalog of real sources but rather to remove all potential sources 
of emission resolved in the shortest wavelength band, an aggres- 
sive source detection level is desirable. Because the confusion 
noise in these maps is already a factor of ~ 3 above the instru- 
ment noise level, lowering the source detection threshold does 
not substantially increase the number of identified sources. The 
250 fim selected sources are then subtracted from the 250 //m 
map to produce the source-subtracted map shown in Eig. [T] In 
order to remove 250jum emission that is correlated with 350 and 
500 //m emission we employ a very conservative approach. The 
first step in this process is to use the 250 //m source candidate 
catalog to determine the possible positions of counterparts in the 
other two channels. The algorithm described below is performed 
independently on both the 350 and 500 yum maps; no direct com- 
parison between the catalogs for the two longer wavelengths is 
performed or necessary. Eor each 250 //m candidate source posi- 
tion, the relevant PRE is fit to the position in the longer wave- 
length channel. In the fit, the position, flux and width of the 
source are allowed to vary. The resulting list of long wavelength 
candidate source counterparts is then compared to the 250 yum 
candidate source catalog. All those fits where the best fit flux is 
different by > 50% at 350yum or > 90% at 500 //m, whose best 
fit width differs from the nominal PRE by > 20%, or whose best 
fitting position is different by a total distance of 2 map pixels or 
more are deemed "unsuccessful". The catalog consisting of the 
successful fits is then subtracted from the long wavelength map 
to produce a source subtracted map at that wavelength. 

Since the angular resolution of the two longer wavelength 
channels is significantly larger than that at 250 fim, it is com- 
mon to have a large fraction of the fits be unsuccessful. As this 



is due to the effects of confusion rather than an actual lack of a 
long wavelength counterpart source, to be conservative we need 
to account for these unsubtracted sources; this process is per- 
formed in two steps. First, we generate a 250 //m map which is 
the sum of the 250/vm source subtracted map (shown in Eig. [T]) 
and all those 250 yum selected sources for which the longer wave- 
length fit was unsuccessful. This corresponds to the 250 //m map 
which would arise if all those sources detected and removed at 
the longer wavelength were also removed from the 250 //m in- 
put map. This map is then convolved up to the longer wave- 
length's resolution using the relevant PRE and rebinned to the 
longer wavelength map's resolution. This map - which is equiv- 
alent to the long wavelength source free 250 jum counterpart map 
- is then pixel-wise fit to the long wavelength map to determine 
the scaling between them. Bright sources tend to make such scal- 
ings difficult in SPIRE data because such sources have individual 
colours which make pixel-wise colour comparisons have large 
variance. However, by construction such bright sources have 
been removed from the comparison so the scaling between the 
colours is well described by a linear model with coefficients of 
0.76 at 350 yum and 0.40 at 500 yum. The scaled 250 yum bright 
source removed map is then subtracted from the equivalent long 
wavelength map. This has the effect of subtracting out both the 
250 jum sources undetected in the longer wavelength map and the 
component of the confused background present in the 250 jum 
map. 

As a further step, the 250 jum sources whose long wavelength 
fits were unsuccessful are binned into a map using their 250 jum 
flux and position, but the PRE of the longer wavelength channel. 
This undetected source model map is then threshold cut; all pix- 
els with values greater than the threshold are deemed possible 
contaminants to the longer wavelength maps and a pixel mask is 
generated based on the positions of the cut sources. The center 
and right most panels of Eig.[T]show the resulting 250/Ym source 
removed and masked maps for the 350 and 500 yum channels. 

The 500 yum map shown in Eig.[T]features an extended source 
coinci dent with the measured center of the SZ effect in this 
cluster iHalverson et alj (l2009l) . If this signal is due to extended, 
roughly azimuthally symmetric emission as expected from the 
SZ effect then taking radial averages of the masked maps will 
enhance the signal to noise ratio of this detection. As the signal 
to noise ratio per map pixel is small, here we assume azimuthal 
symmetry and average the maps in simple radial bins of con- 
stant 18 arcsec width; the results of this are shown in Eig. |2] 
Interestingly, both the 350 and 500 yum data exhibit extended 
emission at the center of the cluster; to constrain this, we fit these 
data to the best fitting isothermal y6 model from Hal verson et al.l 
( l2009h which has a characteristic angular size of a few arcmin- 
utes. The resulting fits are plotted in Eig.|2j Table[T]lists the best 
fitti ng central SZ inc rement brightnesses at 350 and 500 yum. For 
the IHalverson et alj (1200 9) model parameters the probability to 
exceed the measured;t'^ for the fit at 350 fim is 0.2 1 ; for 500 yum it 
is 0.25. The isothermal j6 model is a poor statistical description of 
the radial average at 250 fiva. In addition to the statistical uncer- 
tainties associated with the errors on the data o utput from the fit, 
the uncertainties in the Halverson et alj (l2009l) model contribute 
a further 1.2 MJy sr"' at 350yum and 1.9 MJy sr"' at 500 yum. We 
have not included the contribution of these in the uncertainties 
listed in Table [T] 

In order to compare the SZ effect increments measured here 
with the signal e xpected from previous measurements, we com- 
pile the results of'Andr eani et all d 19991) . IHalverson et al] (l2009h 
and Plagge et al. (2009J. As all of these measurements were per- 
formed with instruments working at 2 mm, they can be aver- 
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Fig. 1. S ource subtrac t ed ma ps at 250,350 and 500 //m; these figures show the same regions defined by the 7'x7'green boxes in 
Fig. 2 of lEgami et alj (l2010h . All three maps are shown on the colour scale shown at right; units are mJy beam"'. The left most 
panel shows the source-subtracted 250 ;um map; the confused sub-mm background is most clearly evident in this band. The center 
and rightmost panels are respectively the 350 and 500 /im source removed and masked maps; their construction is discussed in the 
text. The cross hairs show the best fitting SZ effect centroid from Halverson et al.l (1200 9); the emission visible at 500 //m but not in 
the other maps is consistent with the position, shape and flux from the SZ effect expected at this wavelength. 

Table 1. Best fitting and expected SZ effect parameters at 350 and 500 //m. 



A/o (MJy sr') 
500 jum 350 yum 



A/o(350)/A/o(500) 



SPIRE Best Fit 0.268 ±0.031 0.097 ±0.019 0.36 ±0.12 

Expected with relativistic correction 0.266 ± 0.014 0.060 ± 0.003 0.23 ± 0.05 

Expected without relativistic correction 0.184 ±0.010 0.017 ± 0.001 0.09 ± 0.05 



aged to create a best estimate for the central Comptonization 
parameter using the X-ray determined cluster temperature of 
Te - 13.9keV. Including relativistic corrections, the uncertainty 
weighted average yo is 3.46 ± 0.16 x 10""*. Table [1] lists the ex- 
pected brightness of the SZ distortion in the SPIRE bands in- 
cluding the effect of the actual filter bandpasses; Fig. [3] shows 
the spectral distortion of the CMB brightness for these yo and 
Te values with the SPIRE bandpasses for reference. Good agree- 
ment is found between the expected results given the fiducial yo 
and Te and the SPIRE results at both 350 and 500 yum. 

The relativistic SZ effect corrections required in the pres- 
ence of a relativistic electron population have the largest effect 
at both the peaks of th e SZ effect and the very shortest wave- 
lengths. |Nozawaetal] (2000) has calculated analytic fitting for- 
mulae for these corrections for A > 250 yum; in the calcula- 
tions above we use these formulae to correct the purely ther- 
mal SZ effect spectrum for the presence of relativistic elec- 
trons. To determine the change in the SPIRE result if these fi- 
nite temperature corrections were absent, we also compute the 
SZ spectrum without the relativistic electron correction and cen- 
tral Comptonization parameter which leaves the brightness at 
2 mm constant, yo = 3.16 x 10""^. Fig. [3] shows this pure ther- 
mal SZ effect spectrum and Table [1] lists the corresponding 350 
and 500 jum SPIRE band weighted intensities. 

Table [1] also lists the ratio of the 350 to 500 yum SZ effect 
intensities for the SPIRE measurement presented here, the fi- 
nite temperature SZ effect corrected spectrum, and the pure ther- 
mal SZ spectrum. The SPIRE and relativistic SZ effect corrected 
spectrum are consistent at the lo" level, while the SPIRE mea- 
surement and the purely thermal SZ spectrum are different by 
2.1cr. This is the strongest evidence for the presence of the rela- 



tivistic corrections to date, though not a detection of their effect 
on the SZ spectrum. 



3. Discussion and conclusions 

The possibility of detecting the SZ effect significantly shortward 
of its positive peak is a testament to t he extraordinary ca pabili- 
ties of SPIRE. However, works like Zemc ov et alj (|2007|) show 
that in large galaxy cluster survey samples, significant contami- 
nation to the SZ signal from bright, gravitationally lensed back- 
ground sources is common. Though the properties of this par- 
ticular cluster have not precluded measurement of the SZ effect 
- 1E0657-56 has fewer bright, lensed sources close to its SZ 
effect center than typical clusters which have been observed in 
the sub-mm and is relatively bright and broad in the SZ effect - 
based on this single example it is difficult to determine whether 
gravitational lensing of the sub-mm background will make mea- 
surements similar to this one more challenging in more typical 
clusters at these and other wavelengths. Data from surveys like 
the HLS will allow us to understand whether typical clusters are 
suitable for SZ effect increment detection, and how the lensed 
sub-mm background will affect measurements at other wave- 
lengths. 

The possibility of diffuse emission from cold dust in the 
ICM has been discussed in the past (e.g. IStickel et alll2002h . It 
is expected that, due to sputtering by energetic photons, such 
dust would have a very s hort lifetime in the ICM environment 
dPraine & Salpetenll979h . Based on the radial averages of the 
250jt/m source subtracted maps, where for ICM dust with rea- 
sonable temperatures the brightest thermal emission would oc- 
cur, we find no evidence for this type of diffuse emission in this 
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Fig. 2. Radial averages of the source subtracted and masked 
maps shown in Fig.[T] As the absolute mean of the maps are not 
measured by SPIRE, the points for each band are scattered about 
OMJy sr '. Also plotte d are the best fitt ing isothermal /3 models 
using the parameters in Hal verson et al.l ([2009) in the two SPIRE 
bands in which the SZ effect is expected to be non-zero, 350 and 
500 yum; Table[T]gives the numerical values of the central incre- 
ments brightnesses at both wavelengths. As the 250 /im radial 
bins are not correlated, each point should be compared to the 
(dashed) Hne A/ = OMJy sr ', while in 350 and 500 yum the 
dashed line shows the inferred zero level of the SZ effect. 



cluster Because in any reasonable scenario such dust emission 
would be faint in a z > 0.1 cluster, we expect that targeted 
searches of local clusters will be more successful for this sci- 
ence, though the HLS and similar surveys can provide useful 
constraints. 

Given the statistical uncertainties, the 350 ;um SZ effect in- 
crement measured here is slightly less than Icr larger than would 
be expected from the 2 mm measurements. This may point to 
residual problems with the source subtraction procedure. More 
data will allow tuning of the confused sub-mm background re- 
moval algorithm and checks on whether biases arising from poor 
background removal are endemic in a large cluster sample. 

Given the presence of the SZ effect in clusters, it seems that 
care must be taken with photometry of 500 yum sources within 
an arcminute or so of the cluster center; such sources will be 
positioned on a diff use background so their fluxes will be biased 
by a small amount JRex et al.ll2010l) . 
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Fig. 3. The SZ effect spectrum in 1E0657-56. T he 2 mm uncer- 
tainty weighted average of the mea surements of JAndreani et al.l 
(Il999h . iHalversonetal] (l2009h and lPlaggeet"an (l2009l) (purple 
asterisk) leads to the SZ spectrum shown (solid black curve). 
The best fitting SPIRE measurements at 350 and 500 //m (green 
and red diamonds) and the normalized SPIRE bandpasses (blue, 
green, red solid lines) are shown for reference. The SZ effect 
spectrum which is consistent with the 2 mm measurements but 
excludes the relativistic SZ effect correction is also shown (black 
dotted line). Though both SZ effect curves are well matched 
in the decrement, the finite temperature SZ effect corrections 
change the results by as much as 70% of the expected signal 
in the SPIRE bands. 

As the peak of the FIR emission of dusty sources is 
redshifted to progressively longer wavelengths, high redshift 
sources are expected to have exceptionally red spectral energy 
distributions. The results of this work show that care must be 
taken when searching for sources based on their presence in 
the 500 /im band alone; galaxy clusters whose SZ effects are 
relatively bright and compact could well masquerade as such 
sources. Due to the well known problem of determining coun- 
terparts to sub-mm sources at other wavelengths, such SZ effect 
contamination may not be immediately obvious, particularly as 
cluster fields are crowded and many possible counterparts may 
be present. As a corollary, searching for very red sources in 
confusion limited SPIRE surveys may turn up compact clusters 
based on the presence of strong 500 //m emission. The search 
for such extremely red SPIRE sources in the HLS and other pro- 
grammes is underway now. 
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